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Abstract 


The logistic equation is often considered as a simplification of the 
chemostat (Kooi, Boer, Kooijman (1998)). However, the global quali¬ 
tative properties of the delayed single species chemostat are completely 
known, see e. g. Smith (2011). Such properties still remain open for 
the delayed logistic equation, see e. g. Banheley, Czendes, Kristin, 
and Neymaier (2014) despite that they have been announced a long 
time ago (Wright (1955)). We may therefore ask whether the logis¬ 
tic equation really is a simplification and what information about the 
chemostat actually is contained. We discuss the links between these 
equations and conclude that they are less clear in the delayed case 
than they are in the non-delayed case. We keep the presentation as 
elementary as possible. 

1 Introduction 

In this paper we discuss the relation between the delayed chemostat 


X'(T) 


S\T) = 



X(T -R) - DX(T). (1) 
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and various forms of the delayed logistic equation that are going to be in¬ 
troduced later. Here, S(T) is the substrate concentration at time T, X(T) 
is the concentration of some organism at time T having the substrate S as 
its limiting resource. The parameters C > 0, D > 0, A > 0, B > 0, M > 0, 
R > 0 stand for concentration, dilution rate, search rate for the organism, 
handling time for the organism (cf. Holling (1959)), conversion factor for the 
organism, and time delay, respectively. The time delay could be included 
in many ways in the chemostat model, but we think that the above way is 
the most natural way to include it; in this way it takes into account that 
the organism cannot convert consumption of the substrate into own biomass 
instantaneously. Note the factor e~ DR , it corresponds to the fact that a part 
of the individuals that consume nutrient is washed out before they are able 
to reproduce. 

The model above is a phenomenological model for a number of ecological 
phenomena, cf. Smith and Waltman (1995). For instance, a lake ecology may 
have rivers transferring a limiting resource to the lake and rivers diluting this 
resource. 

We start our analysis by re-parameterizing our model into a dimensionless 
form by the changes t = DT , r = DR, s = S/C, x = X/C, a = AC/D, 
b = ABC, and m = M. For simplicity, we introduce the function / defined 


by 


f(s) 


as 

1 + bs' 


too. We note that / is strictly increasing with /(0) = 0. We arrive in the 
dimensionless form 


s(t) = 1 - s(t) - f(s(t))x(t), 

x(t) = me~ r f(s(t — r))x(t — r) — x(t) (2) 

for which we impose the initial conditions s(t) > 0, — r < t < 0 and x(t) > 0, 
—r < t < 0, x(0) > 0. A number of important results regarding the model 
(]2|) has been proved and can be found in Smith (2011). They are that the 
solution space is C([—r, 0], R 2 ), solutions remain non-negative (Theorem 3.4), 
and bounded (Section 5.6). Either one or two equilibria, the washout state 
and the survival state exist, the washout state is globally attracting if it 
is the only equilibrium and so is the survival state it it exists (Theorem 
5.16). The Erst proof of the global stability result can be found in Ellermeyer 
(1994). Thus, all global qualitative dynamical properties of the single species 
chemostat model are known. 
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2 Methods 


We briefly go through the main methods that are going to be used in this 
paper. Most of this material can be found in the monographs of Hale and 
Lunel (1993), Smith (1995), and Smith (2011). It turns out that ((2]) can be 
reduced to a scalar equation, so we work throughout the paper with scalar 
autonomous delay differential equations of the form 


x(t) = f(x(t),x(t - r)) 


(3) 


satisfying an initial condition x(t) = (j>{t), —r < t < 0. The solution of 
(J3J) starting at the initial condition (j) G C([— r, 0],R) is denoted by x(t,cj)). 
We denote the right hand side by f(x,y ) when we want to refer to the 
different variables or take partial derivatives. Existence and uniqueness of 
the solutions of such equations are granted by the method of steps if / and f x 
are continuous on R 2 and (j) is continuous on — r < t < 0 since the equation 
becomes an ODE after substituting the continuous initial condition into the 
equation. Hence, standard ODE theory applies and grants existence and 
uniqueness of solutions of equations of form ((3]) . It follows that the natural 
state space of the equation ([3]) is the continuous functions on the interval 


r < t < 0, or C([—r, 0], R). 


The equations that we are going to consider arise in population dynam¬ 
ics and we are going to use a theorem granting that non-negative solutions 
remain non-negative. Such theorems require that solutions exist and are 
unique, but additional conditions are required. This additional condition is 
given by f(0,y) > 0, Vt G R,Vy G R+ and we note that this condition 
can be expressed in terms of the functions involved in the differential equa¬ 
tion directly. Therefore, direct substitution in the involved expressions can 
be used to check these things. Also this theorem is directly related to ge¬ 
ometric ODE theory, it basically formulates conditions for when the vector 
field is either directed into the positive cone or along the boundaries of it. 
No problems arise when it is directed into the positive cone and in the latter 
case uniqueness of solutions grants that no solutions escape from the positive 
cone. 

Studies of nonlinear ordinary differential equations usually start from 
analysis of equilibria. All equilibria can be found by neglecting delays and 
therefore they coincide with the equilibria of the corresponding ODE 


x(t) = f(x(t),x(t)). 


(4) 
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We claim that we have found one equilibrium by solving f(x,x) = 0, say 
x(t) = x. The linearization of ((2D around the equilibrium x is then given by 

i(t) = fx(x, x)£(t) + f y (x, x)£(t - r) = a£(t) + 6£(t - r). (5) 

We shall as usual look for exponential solutions around the zero solution of 
this equation and these are characterized by the characteristic equation 

A = a + be~ Xr . 

The solutions of © classifies the solutions of © into four cases according to 
the following theorem (Theorem 4.7 in Smith (2011)): 

Theorem 1. The following hold for ()5lh 

(A) If a + b > 0, then £ = 0 is unstable 

(B) If a + b < 0, b > a, then £ = 0 is asymptotically stable regardless of the 

magnitude of the delay. 

(C) If a + b < 0, b < a, then there exist threshold value for the delay r* 
(which can be explicitly computed) such that £ = 0 is asymptotically 
stable for 0 < r < r* and unstable for r > r*. 

(D) If a + b = 0, then A = 0 is a root of ((2j) . 

In case (D) the linearization ((5j) does not provide the qualitative infor¬ 
mation requested regarding solutions of the original nonlinear equation ((3j) 
in the vicinity of its equilibrium x. However, center manifold theory (Guck- 
heimer and Holmes (1983)) exist in the delay differential equation case, too, 
see Diekmann, van Gils, Verduyn Lunel, and Walther (1995). 

In ODE theory there exists a standard method for constructing Lyapunov 
functions for linear equations when all characteristic roots have negative real 
part, see Jordan and Smith (1990). There are extensions of this method to 
delay equations, see Hale and Lunel (1993), but in most cases explicit con¬ 
structions can be very complicated already in simple cases and most results 
exist just in the case (B) above. To get Lyapunov functions that are valid for 
the nonlinear equations like (|3|) . such results still usually need to be modified. 

The method of Lyapunov functions and LaSalle’s (1960) invariance prin¬ 
ciple is always the most desirable when proving global stability. However, 
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more straightforward methods might exist. We start by mentioning the fol¬ 
lowing simple comparison theorem (Smith (2011), Theorem 3.6). Essentially, 
it states that if a solution start in front of (behind) another and grows faster 
(slower), it will stay in front of (behind) the other. 

Theorem 2. Consider (ED and suppose that f and f x are continuous on R 3 
and that the initial condition </> is continuous on [— r, 0] . Let x(t) be a solution 
of fl3D with x(t) = (j>{t) for [— r, 0] on some interval [0, s], 3 > 0. Let x(t) be 
continuous on [— r, 0] U [0, s], differentiable on [0, s] and satisfy 

x(t) > (<) f(x(t),x(t-r)), t G [0, 3], 
x(t) > (<) 4>(t),—r<t< 0. 

Then x(t) > (<)x(t). 

We continue formulating conditions that preserve the order of two solu¬ 
tions as follows. It is Theorem 5.9 in Smith (2011). 

Theorem 3. Consider (ED- Assume that f : R 2 —R, f x , f y are continuous, 
f y (x,y) > 0, and that x\(t) and x- 2 (t) are two solutions of (J3D defined on 
[— r, s] ; for some 3. If x\{t) < x^{t) for t e [— r, 0], then x\(t) < X 2 (t) for 
t E [-r, 3]. 

This theorem is a direct consequence of the previous theorem and makes 
it possible to translate many of the properties of the equilibria of the scalar 
ordinary differential equation (ED to (ED- The following theorem contains the 
parts of Corollary 5.11 in Smith (2011) that we shall need. 

Theorem 4. Let x G (a, b) be an equilibrium of (ED, and hence of (ED- Sup¬ 
pose that the sign-condition 

(x — x)f(x,x)< 0, xE[a,b \, x^x 

holds. If 0(s) G [a, b], s G [— r, 0], then 

lim = x. 

£—>•00 
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3 A delayed hyperbolic model 

A relation between the single species chemostat model and the logistic model 
exists, see e. g. Kooi, Boer, Kooijman (1998) and this relation can be made 
visible also in the delayed case. More precisely, the logistic equation can be 
considered as an approximation of the mass-balance equations used in the 
chemostat. Let us consider the functional 


V(t ) = x(t) + me r s(t — r) — me 


( 6 ) 


Along the solutions of (j2J), we identify that equation for the total derivative is 
the ordinary differential equation V(t) = —V(t) and thus, lim^oo V(t) = 0. 
We ask what models we could obtain by considering the delayed chemostat 
in the hyperplane V = 0 in C([— r, 0], R 2 ). It is not evident that such a 
procedure will preserve the dynamics even in the finite dimensional case 
(Tliicme (1992)) but we try this substitution and get 


x{t) 


me r f 


me 


x(t) 


me 


x{t — r) — x(t). 


(7) 


This is a single hyperbolic delayed equation. Our program is for the moment 
to verify that this equation has the same qualitative dynamical properties as 
the single species chemostat model. Indeed, we prove the following sequence 
of results. 


Lemma 1. Solutions o/([7J) remain non-negative and bounded above by x(t) < 
me~ r . 


Proof. We start be proving that the solutions remain non-negative. This 
follows from the fact that the estimate 


me r f(l)x(t — r) > 0 

holds for all positive x(t — r). Theorem 3.4 in Smith (2011) gives the result. 
Comparison with 

x(t) < —x(t) 

for x(t) > me~ r shows that large positive solutions decay at least exponen¬ 
tially since f(s) < 0 for s < 0. It follows that positive solutions are bounded 
above by x(t) < me~ r . □ 
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Lemma 2. The equation (J7J) has the washout state x = 0 as its unique non¬ 
negative equilibrium if /(l) < —. This equilibrium is then locally stable. If 
/(1) > z~, then ([7|) has two non-negative equilibria, the washout state that is 
locally unstable and the survival state that is locally stable. 

Proof. By neglecting delays, we get that the washout state is an equilibrium. 
For the survival state we solve 

/ me~ r — x{t) \ e r 
\ me~ r ) m 

and conclude that this equation has a unique positive solution x{t) — x if 
and only if /(1) > —. In order to investigate the local stability of these 
solutions, we compute the characteristic equation of a generic equilibrium of 
the equation (JTj). It takes the form 

A + 1 + xf'(s) — e~ Xr me~ r f(s) = 0 

that turns out to be the second factor of the characteristic equation of the 
delayed single species chemostat (j2J), see Smith (2011, p55). The factorization 
of this equation given there is therefore, not a coincidence. It is closely related 
to motion in the hyperplane V = 0 defined by (jSJ) and motion towards it. 
Hence, we have proved that the local stability properties of the delayed single 
species chemostat (J2]) and the hyperbolic model ([7j) are the same. □ 

We go on formulating the global properties of (J7J) and remark that we use 
the monotone dynamics of it in order to arrive in the conclusions. This is a 
simplification in comparison to proving the global stability properties of the 
chemostat equations (j2|) directly, since the first equation does not satisfy the 
quasimonotone condition in Smith (2005). 

Theorem 5. The following statements hold for the solutions of 0. 

(i) ///(1) < then the washout equilibrium is the unique non-negative 
equilibrium of (J7]) and it attracts all non-negative solutions. 

(ii) ///(1) > —, then the survival equilibrium exists for (13) and it attracts 
all positive solutions. 

Proof. Consider (171). By Lemma [TJ we have 0 < x(t) < me~ r . With aid of 
the differential inequality 

x{t) < me~ r f(T)x(t — 1) — x(t) 
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we notice that no solutions can grow faster than the solutions of the linear 
delay differential equation 


z(t) = me~ r f{l)z{t — 1) — z{t) (8) 

which has the characteristic equation A = —1 + me~ r f(l)e~ x . Now 

me"7(l) < 1 (9) 

implies that the zero solution of ([S]) is asymptotically stable and falls into 
category (B) of Theorem [0 Therefore, all solutions of (J7J) converge to zero 
if (jUJ) since they are bounded from above by solutions that decay to the zero 
solution. 

Now assume me~ r 1) > 1 and remember that x(0) > 0. The differential 
inequality x(t) > —x(t) gives x(t) > 0. We can therefore, work with initial 
conditions satisfying x(t) > A, for some A > 0 and — r < t < 0. Consider (J7]) . 
It possesses monotone dynamics (or order-preserving dynamics, see Smith 
(1995)) since 0 < x{t) < me~ r gives 


me r f 


me r — x(t ) 
me~ r 


> 0 


so that Theorem [3] can be applied. Consider the non-negative equilibria of 
(J?P that are 0 and x when me~ r f( 1) > 1. Both equilibria are shared with 


the ODE 


x(t) = ( me r f 


me 


x(t) 


me" 


1 x(t). 


For the above ODE, the survival state is globally asymptotically stable on 
x{t) G (0,me _r ). Theorem |4] gives that x(t) = x attracts all solutions of (171) 
with initial conditions satisfying A < x(t) < B with A > 0 and B < me~ r . 
In the beginning of this proof and in the proof of Lemma [T] we have justified 
that we can work with such initial conditions provided x(0) >0 and x(t) > 0 
for — r < t < 0. □ 


Thus, the delayed hyperbolic model (J7J) preserves the dynamics of the 
delayed single species chemostat model ((2j). We conclude that 6 = 0 gives 
f(s) = as and therefore the delayed logistic model 


x(t) = ame r x{t — r) — x(t) — ax(t)x(t — r ). 


( 10 ) 




shares the same dynamical properties, too. In particular, it has monotone 
dynamics when 0 < x(t) < me~ r . If we let the delay parameter tend to zero, 
this equation becomes the ordinary differential equation 

x(t) = (am — 1 )x(t) — ax 2 (t ) = (am — 1 )x(t) I 1 — 


We can identify the growth rate and carrying capacity parameters as am — 1 
and (am — 1 )/a, respectively. The washout equilibrium is still given by 
x(t) = 0 and the survival equilibrium is located at the carrying capacity 
x(t) = (am — 1 )/a. 

4 The logistic equation 

Delayed logistic equations are not always derived in a way that preserves the 
dynamical properties of the underlying chemostat equations. Consider, for 
instance, Hutchinson’s (1948) delayed logistic model 

x(t) = (am - l)x(i) (l - X ^ m _P ) ■ (11) 


We remind the reader that the resource dynamics was not modeled explic¬ 
itly in Hutchinson’s paper and that the problems mentioned by Kooi, Boer, 
Kooijman (1998) exist. This equation can be made comparable to (TTCTT) by 
arranging it as 



x(t) = am ■ x(t ) — x(t) — ax(t)x(t — r). 

We note that the above equation has preserved the equilibria at x(t) = 0 and 
x(t) = (am — 1 )/a. However, it does not have monotone (order-preserving) 
solutions, since 

— ax(t ) < 0. 

We start by linearizing around a generic equilibrium x and get the charac¬ 
teristic equation 

A = (am — 1 — ax) — e~ rX ax. 

For x(t ) = 0, we have A = am — 1 and for x(t) = (am — 1 )/a, we have 

A = —e~ rX (am — 1) 
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which according to Theorem [Q belong to the category (C) above: There exists 
a Hopf-bifurcation value for the magnitude of the delay, so that the equilib¬ 
rium is stable if the delay is below that threshold and unstable otherwise. 
We have now several indications of the fact that the dynamics of (ITT]) is 
considerably more complicated than the dynamics of (TTQj) and this is not the 
end of the story. 

In order to investigate (fTTD further, we introduce r = (am — l)t, p = 
(am — l)r and 

£ = — - 1 

s am —1 

a 

We get 

^ = i+m) ( 12 ) 

which is the famous Wright’s (1955) equation. The following is known about 
(H3- First, its solutions are bounded (Proposition 5.13 in Smith (2011)). 
For 0 < p < 1.5706, every orbit tends to the survival state (Wright (1955) 
and Banheley, Czendes, Kristin, and Neymaier (2014)). Note that 1.5706 < 
| « 1.5707963268 and that Wright’s (1955) conjecture claims this to be 
true for all values of the delay in 0 < p < |. This problem has been open 
for almost 60 years now and other open problems connected with Wright’s 
equation exist, too, cf Banheley, Czendes, Kristin, and Neymaier (2014)) and 
references therein. We can therefore ask whether it is really a simplification 
to use logistic equations instead of the chemostat when modeling the growth 
rate of the lowest trophic level. The single-species delayed chemostat has a 
complete qualitative analysis. There exists at least one periodic solution for 
(112]) when p > f. 

5 Summary 

In this paper we have studied a delayed model for the single species chemo¬ 
stat. The global dynamics of this model has been completely known since 
Ellerrneyer (1994) and new variants of such an analysis has been given after 
this (Smith(2011)). All solutions tend to either the washout or the survival 
state and the local stability properties of those equilibria fully determines 
the global properties of the model. 

We make a new interpretation of these results here. The results can be 
divided into motion towards a hyperplane in C([— r, 0], R 2 ) and motion in the 
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hyperplane governed by a scalar monotone delay differential equation. We 
name this equation, the hyperbolic model. Such a separation is not always 
valid even in the finite dimensional case (Thieme (1992)). We verify that 
no qualitative differences occur between the chemostat equations and the 
hyperbolic model. 

A limiting case of this hyperbolic model correspond to a delayed logistic 
model. We prove that the delayed logistic growth rate model encountered 
in this way does not correspond to the model that usually is referred to as 
the delayed logistic model in the literature (Hutchinson (1948)). This model 
does not satisfy any monotonicity conditions (Smith (1995)) and is equipped 
with a number of problems regarding its long-run dynamical behavior that 
has been open for a long time, see Banheley, Czendes, Kristin, and Neymaier 
(2014). 
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